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ABSTRACT 

Massive stars generally end their lives as neutron stars (NSs) or black holes (BHs), with NS formation 
typically occurring at the low mass end and collapse to a BH more likely at the high mass end. In an 
intermediate regime, with a mass range that depends on the uncertain details of rotation and mass loss 
during the star's life, a NS is initially formed which then experiences fallback accretion and collapse to 
a BH. The electromagnetic consequence of such an event is not clear. Depending on the progenitor's 
structure, possibilities range from a long gamma-ray burst to a Type II supernova (that may or 
may not be jet-powered) to a collapse with a weak electromagnetic signature. Gravitational waves 
(GWs) provide the exciting opportunity to peer through the envelope of a dying massive star and 
directly probe what is occurring inside. We explore whether fallback onto young NSs can be detected 
by ground-based interferometers. When the incoming material has sufficient angular momentum to 
form a disk, the accretion spins up the NS sufficiently to produce non-axisymmetric instabilities and 
gravitational radiation at frequencies of ~ 700 — 2400 Hz for ~ 30 — 3000 s until collapse to a BH 
occurs. Using a realistic excess cross-power search algorithm, we show that such events are detectable 
by Advanced LIGO out to w 17Mpc. From the rate of nearby core-collapse supernovae, we estimate 
that there will be ~ 1 — 2 events each year that are worth checking for fallback GWs. The observation 
of these unique GW signatures coincident with electromagnetic detections would identify the transient 
events that are associated with this channel of BH formation, while providing information about the 
protoneutron star progenitor. 

Subject headings: black hole physics — gravitational waves — stars: neutron — supernova: general 



1. INTRODUCTION 

Determining the fate of zero-age main-sequence 
(ZAMS) stars with large masses is a long-standing prob- 
lem. In general, it is expected to depend in a compli- 
cated way on mass loss and rotation during the star's 
life. These in turn are related to details such as the mag- 
netic field, metallicity, and binarity. Even with all these 
uncertainties, theoretical efforts indicate a rough general 
picture. A core-collapse supernova that successfully un- 
binds its stellar mantle leaves a neutron star (NS) behind. 
In cases when this does not happen, a stellar-mass black 
hole (BH) is instead expected, but this can occur in a 
nu mber of different ways; ( for a more detailed discussion, 
see IQ'Connor fc Qtti 1201 lh . For example, if there is a 
nuclear phase transition during protoneutron star (PNS) 
cooling, or if cooling reduces pressure support in a hyper- 
massive PNS, a BH results. In another scenario, if the 
supernova mechanism fails to revive the accretion shock, 
continued accretion pushes the PNS over its maximum 
mass, creati ng a BH with likely lit tle or no electromag- 
netic signal (jKochanek et al.1 120081 ). Finally, if the core- 
collapse supernova is successful, but perhaps weak, then 
the young NS will be subject to fallback accretion rates 
of M ~ 10~ 4 - 10~ 2 Mp) s" 1 over the next ~ 30 - 1000 s 
(jMacFadven et all [20011 IZhang et al.ll200l . This addi- 
tional material pushes the NS past its maximum mass, 
again resulting in a BH. 

In the present work, we focus on this latter fallback 
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mechanism for creating BHs. Si nce the i d ea of fallback 
accretion was first discussed by iColgatej (|1971l) , it has 
been an important area of focus for theoretical stud- 
ies (e.g.. lChevalier|[l989[ iWooslev fc Weaver| [l995t iFTverl 
[l99alZhang et al.ll2008HUgliano et al.ll2012h . There have 
also been a wide range of predictions for the type of 
events associated with fallback accretion leading to BHs. 
When rotation is included, the newly formed BH con- 
tinues to accrete and may produce a jet. Depending on 
the mass of the envelope at the end of the star's life, 
this may result in a jet -powered Type II supern ovae or a 
long gamma-ray burst (jMacFad ven et al.ll2001l ). In cases 
where the jet is pointed away from the observer or jet 
formatio n does not occur, a dim supernova could re sult 
instead (jFrver et all [20071 120091 iMoriva et al.l [2010h . If 
a quickly-spinning, strongly-magnetized NS is present, 
it may fling away the in falling material an d produce a 
"propeller nova" instead (|Piro fc Qttl(20TTh . 

Observationally, it would be helpful to determine the 
main-sequence mass range that leads to BHs created 
via fallback, so tha t it can be compar ed with theoret- 
ical expectations. iHeger et al.1 (j2003| ) argue that for 
sub-solar metallicities, this occurs for a ZAMS mass 
of ~ 25 — 40M Q , but recent work paint s a more 
complicated picture (jO'Connor fc Ottl 120 111) . It has 
been well-established that most long gamma-ray bursts 
(which may follow fallback accre tion) are associated with 
broad-lines Type Ic supernovae (IWooslev fc Bloomll200a 



iHforth fc Bloomll2011l : lModiazll201lt ). but such events are 
too distant to directly identify the exploding stars. Pro- 
genitor stars associated with standard Type II-P super- 
novae via pre-s upernova imaging generally have masses 
17 - 20M o (jSmartt et al.ll2009D . which is lower than 
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the maxi mum mass expecte d for fiducial core-collapse 
scenarios (|Heeer et al.ll200l . This might suggest that 
other types of supernovae, and maybe even BH for- 
mation, begin occurring in a lower mass range. Such 
a prospect is maybe not surprising given the histori- 
cal difficulties in robustly producing supernovae via the 
neutrino mechanism in theoretical models ( although see 
iMueller et aTll2012[ ). On the other hand, iSmith et all 
( 20111) argue the mass range of Type II- P progenitors 
are consistent with a substantial fraction (~ 25 — 35%) 
of supernovae events being produced in binary systems, 
so the situation remains unresolved. 

Gravitational waves (GWs) provide an independent 
probe to determine what processes are occurring deep 
within these mass ive stars (for example, see the stud- 
ies and r eviews bvlFrver et alJl2002HDimmelmeier et al.1 
I2002allbl : lOtt et all 120041 IQttl 1200ft IKotakd 120 lift . The 
general picture we explore is as follows. Assuming that 
the fallback material forms a disk before reaching the 
NS, the NS accretes sufficient angular momentum that 
its spin parameter /3 = T/\W\ reaches a critical value j3 c . 
Here T is the rotational energy and W is the gravita- 
tional binding energy. Above /3 C , non-axisymmetric in- 
stabilities occur and GWs are radiated (the exact value 
of /3 C depends on a number of factors, which are discussed 
in detail in H2.3|) . Since the star is quickly torqucd down 
when j3 > /3 C , and quickly spun up by accretion when 
j3 < /3 C , the NS is forced into a state of marginal insta- 
bility with p sa j3 c while it continues to gain mass. The 
result is ~ 30-3000 s of high frequency (~ 700 -2400 Hz) 
GW production until the NS becomes sufficiently mas- 
sive to collapse to a BH. As we describe in more detail 
below, the detection of such a GW signal would be strong 
evidence that fallback accretion is occurring. 

In SJ21 we summarize the fallback model that is em- 
ployed as well as our treatment of the non-axisymmetric 
shape of a quickly spinning NS. In JJJ we provide a de- 
tailed discussion of the detection te chniques using a ex- 
cess cross-power search algorithm (|Thrane "eTaTT[20Tll) 
to recover an example waveform. Given the inherent 
uncertainties in the messy process of fallback accretion, 
we consider this a much more representative assessment 
of the detectability than a naively optimistic matched- 
filtering approach, and we estimate that aLIGO can see 
such events out to « 17Mpc. In fj4] we discuss theoreti- 
cal and observational constraints on the types of electro- 
magnetic transients expected to be associated with such 
events. We conclude in fj5] with a summary of our re- 
sults, a discussion of future work, and some speculations 
about what astrophysics can be learned from these GW 
detections. 

2. THE FALLBACK AND SPINNING NS MODELS 

We begin by summarizing the main features of our 
models for fallback accretion and our treatment of spin- 
ning NSs. These semi-analytic models are clearly ideal- 
ized, but they allow us to survey parameter space for this 
initial work. In future investigations, we plan to employ 
more realistic models based on numerically calculated 
progenitors. 

2.1. Fallback Accretion 

We conside r fallback accretion onto a newl y born NS as 
discussed by iMacFadven et all (|200lD and IZhang et all 



(|2008l) . In particular, we focus on a progenitor with a 
25M ZAMS mass with relatively ineff icient semiconvec- 
tive m ixing (referred to as model A in IMacFadve n et all 
2001; . This star experiences mass loss during its lifetime 
and i s reduced to 14.6-M^ by the time core collapse oc- 
curs. IMacFadven et all (|2001l ) input by hand a range of 
explosion energies from 2.6 x 10 50 — 1.3 x 10 51 erg and 
find a range of fallback masses from 0.24 — 3.7 Mq, where 
the two are inversely related. In all cases the remainder 
of the star that does not fall back is unbound. For ener- 
gies greater than this, the entire envelope is ejected and 
a NS remnant is left. 

To replicate the main features expected for fallback, we 
approxi mate the accretion rate as two power-laws (as was 
done in lPiro &; Otdl20Tlh . At early times it is relatively 
flat and scales as 

Mcarly = ^IQ-H^Mq B~\ (1) 

where 77 « 0.1 — 10 is a factor that accounts for different 
explosion energies (a smaller 77 corresponds to a larger 
explosion energy), and t is measured in seconds. This 
scaling of 77 by two orders of magnitude corresponds to 
a change of the explosion energy from 2.6 x 10 50 — 1.2 x 
10 51 erg, which is merely a factor of five difference. This 
demonstrates just how sensitive the accretion rate is to 
the explosion energy. The late time accretion is roughly 
independent of the explosion energy and just depends on 
the mass of the progenitor at the onset of collapse. It is 
set to be 

Af latc = 5Ot" 5 / 3 Af s- 1 . (2) 
Interpolating these two expressions, we use 

M=(M~X + M^ e y 1 (3) 

for the accretion rate at any give time 

For a NS with initial mass Mo, the time-dependent 
mass is 

M{t) = M + f Mdt. (4) 

JO 

Note that M(t) corresponds to the total baryonic mass o f 
the NS, but as discussed in iLattimer fc Prakashl (|2001f) . 
a non-negligible fraction of this mass becomes binding 
energy and is radiated away in the form of neutrinos. 
For a baryonic mass Mb ary on, the gravitational mass of a 
remnant with radius i? gra v is 

M grav = M baiyon (l + l G ^ h ^T ) 1 . (5) 

\ -Kg rav C J 

Depending on the mass and radius, this can amount to 
a « 5 — 30% correction. Given that we are just roughly 
approximating the accretion rate and that the conversion 
from a baryonic mass to a gravitational mass will occur 
in a time-dependent manner, we ignore this complication 
for this initial study. 

2.2. Spinning NS Model 

Besides increasing the NS mass, accretion also causes 
the NS to gain angular momentum. We assume that the 
infalling material roughly circularizes before reaching the 
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NS surface. For this to occur, it must have specific an- 
gular momentum j > (GM i? ) 1/2 ps 2 x 10 16 cm 2 s -1 , 
where Rq ~ 20 km is the nonrotating radius of the young 
NS. For material initially at a radius of 10 10 cm, this cor- 
responds to a rotational velocity of merely ~ 20 km s _1 , 
so it appears quite likely that the NS accretes from a 
disk. The torque the NS experiences is therefore 



A, 



M(GMR ( 



a/2 



(6) 



where R e is the NS radius at the equator. We use R e to 
differentiate it from the nonrotating radius Rq since in 
general R e > Rq for a rotating body. 

Equation (|6]) uses the accretion rate directly from fall- 
back of the envelope, but more exactly this accretion rate 
should reflect the rate that mass is transferred through 
the disk. To explore whether this leads to a quantitative 
change of the accr etion rate, we built one-zone, a-disk 
models (similar to lMetzger et al.ll2008D using the angu- 
lar mo mentum profiles of the m assive, rotating progeni- 
tors of [Wooiiey^jlilil] ([2006|) . Our general finding was 
that the disk reaches nearly steady state, where the ac- 
cretion rate onto the star differs from the infall rate by 
no more than a factor of ~ 5 (and this scales with the 
a- viscosity, with a larger a resulting in higher accretion 
rates). We therefore consider the mediation of the disk 
to be degenerate with rj and use the direct infall rates 
as described above. In future studies using infall from 
numerical simulations, we plan to include these effects of 
a disk. 

The rotation rate of the NS is measured in terms of the 
spin parameter (3 = T/\W\. In the absence of an energy 
loss mechanism, the equilibrium shape of the spinning NS 
is characterized as an axisymmetric, Maclaurin spheroid. 
The relat ion between and t he ellipticity e for these 
figures is (|Chandrasekharlll969ft 
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eil-e 2 ) 1 / 2 



(7) 



where e 2 = 1 — (R z /R e ) 2 and R z is the vertical radius 
oriented along the spin axis. Although this result is for 
an incompressible fluid, it is also rough ly valid for com - 
pressible configurations for our use here ((Lai et al.lll993h . 
since we never consider anything that is more than very 
mildly triaxial. 
The spin of the NS is given by 



2nGp 



(1-e 2 ) 1 / 2 



(3 - 2e 2 ) sin -1 e 



3(1 -e 2 



where p = 3M/AttRq is the average density and q n = (1 — 
n/5)n n with n as the polytropic in dex and k„ as a con- 
stant of order unity (see Table 1 in lLai et al.| [T993?) . The 
dimensionless solutions for Maclaurin spheroids given by 
equations and ((HJ arc summarized in Figure [T] Also 
shown in the bottom panel as a dashed line is the fre- 
quency dependence of the Jacobi ellipsoids, which are 
mentioned later. These must be found numer i cally, and 
the plotted solutions are from iChand rasckharl (|1962f) . 

The mean radius of the rotating star R = (RzR 2 ) 1 ^ 
is given by 



R — Rq 



sin 



:(l-e 2 )V 6 (l-/3) 



-Thj (3 — n) 



(9) 
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Fig. 1. — General solutions for a Maclaurin spheroid given by 
equations 10 and JSJl. The vertical dotted line indicates the lo- 
cation of /9 = /3 aec = 0.14. The dashed line in the bottom panel 
shows the ellipsoidal Jacobi solutions that branch off at /3 Bec . 
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This shows how the equatorial radius may be dramat- 
ically increased when the eccentricity is high, which in 
turn increases the specific angular momentum of accreted 
material. 

In Figure [2] we plot sequences of Maclaurin spheroids 
for a NS with M = 1.3M or M = 2.5M© (as labeled in 
the upper panel) and Rq = 20 km, comparing polytropic 
indices of n = 0.5 (dashed lines) with n = 1 (solid lines). 
This brackets the range of reasonable values for NSs. In 
the top panel we plot the spin frequency as a function of 
f3. As the mass increases from 1.3M© to 2.5M Q , the cor- 
responding spin frequency increases at fixed (3. Similarly, 
the more centrally concentrated n = 1 polytrope gener- 
ally exhibits a higher spin frequency than the n = 0.5 
polytrope. In the bottom panel we plot both the average 
radius R and the equatorial radius R e . The equatorial 
radius is strictly increasing as a function of /3, demon- 
strating how the shape is becoming increasingly like a 
flattened pancake. Changes in the average radius are 
more modest. 

One issue that deserves mention is that for n = 1, solid- 
body rotating models have equatorial velocities that ex- 
ceed Keplerian at merely (3 ps 0.08. This means that 
such models are not self-consistent, but they can be sta- 
bilized against equatori al mass loss with a sm all amount 
of differential rotation (lLai fc Shapiro! fl995l ). For now 
we ignore this detail, but also consider n — 0.5 poly- 
tropes that are at least stable against mass shedding up 
to 13 ps 0.15. 

2.3. Instability and GW Production 

Although these equilibrium Maclaurin spheroids can 
be found for all /? plotted in Figure [H such configura- 
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Fig. 2. — Spin frequency fi, average radius R, and equatorial 
radius i? e as a function of the spin parameter 0, comparing n = 0.5 
(dashed line) and n = 1 (solid line). The NS has either M = 1.3M Q 
or M = 2.5Mq (as labeled in the upper panel), but in either case 
Ro = 20 km. The critical /3 BEC = 0.14 for secular instability is 
plotted as a vertical, dotted line. 

tions are not in general stable, and liable to break their 
symmetry and transition to a nonaxisymmetric shape. 
When /? exceeds a some critical (3 C and instability sets 
in, the resulting triaxial shape causes the production of 
GWs. This in turn limits /3 from any further growth. On 
the other hand, if /3 is forced below (3 C , GW production 
ceases, and the NS spins up again. Due to these compet- 
ing effects, we expect the NS spin parameter to reach a 
saturated state where j3 « f3 c . The NS then continues to 
accrete as it produces GWs at a rate that maintain this 
balance. 

There are a number of potential values for (3 C that could 
be considered. For example, it has been well-established 
that a dynamical bar-mode instability sets in when j3 > 
Pdyn = 0.27. This lea ds to mass shedding and spindown 
back to a stable state (jShibata et al.l l200tf). At lower val- 
ues of j3, but still greater than (3 SCC = 0.14, the possible 
solution for the spinning NS can also be triaxial and given 
by a Jacobi ellipsoid (as plotted in Figure [T]) or Dedekind 
ellipsoid (not plotted since formally these figures are not 
rotating). These represent lower-energy configurations, 
but the NS can only trigger these so-called secular insta- 
bilities and transition from a spheroidal Maclaurin solu- 
tion to these solutions if acted upon by some sort of dis- 
sipative process. Dynamical shear instabilities may also 
operate for (3 > 0.01 if differential rotati o n is present 
(iCentrella et all 120011; iShibata et all 120031: iWatts et al.l 
l2005HOu fc Tohlindl2006HCorvmo et al.ll2010[) . Since the 
conditions required for these low j3 instabilities and their 
associated GW signature are more complicated, we focus 
on the secular instability at /3 sec for the majority of this 
work. 

Potentially destabilizing mechanisms that may trigger 
instability for j3 > /3 SCC , and that have been well-studied 
in the literature, arc viscosity and gravitational radia- 



tion r eaction (iChandrasekharl 119701: iFriedman fc Schutzl 
[19781 lLaHf2001t IGaertig fc Kokkotasll201lD . A rough es- 
timate for the growth time due to dcstabilization from 
gravitational radiation is 



/ gw 



2 x IQ- 5 M^R 4 20 (f3 



/^sec) 



-5 



(ii) 



which is roughly independent of n (jLai &: Shapirol[l995( ) . 
This destabilizes the Dcdekind-mode, which corresponds 
to a highly differentially rotating, but stationary figure 
with a rotation pattern that gives rise to a bar-mode like 
oscillation. In comparison, typical timescales for accre- 
tion arc r acc M/M <~ 10 2 — 10 4 s. For gravitational 
radiation reaction to be effective, it must therefore act in 
a timescale r gw < r acc , which implies (3 > 0.16 — 0.18. Al- 
though not far above f3 scc = 0.14, it is possible that other 
destabilizing processes can occur first as the NS spin is 
increasing. For example, the accretion itself or viscous 
processes near the surface of the NS may be extremely 
dissipativc. For these reasons we consider it reasonable 
to assume that /3 C w /3 SCC , but whether there are small 
variation in f3 c for the particular case of fallback accretion 
deserves further study in the future. 

The next important question is what kind of rotating 
figure is present during this saturated state. In general, 
there are a wide range of spheroidal and ellipsoidal fig- 
ures available, but for the moment we concentrate on the 
Dedekind and Jacobi ellipsoids. Although these belong 
to the much larger class of Riemann-S ellipsoids, the gen- 
eral features of our argument still apply. The Dedekind- 
mode is destabilized by GWs because of its stationary 
properties. This allows it to have an equilibrium state 
without GW production, and at the same time it exhibits 
high internal differential rotation to conserve its circula- 
tion. In a sense, the stable properties of the Dcdekind- 
mode are somewhat unique to the effect of GWs, which 
only care about the time- varying quadrupolc moment . In 
contrast the Jacobi-modc is destabilized for any viscosity. 
Furthermore, its rigidly rotating interior (with perhaps 
a small amount of differential rotation to prevent mass 
shedding) seems more representative of the interior of a 
NS, which is subject to a wide range of potential vis- 
cosity mechanisms as well as magnetic fields that resist 
any differential rotation. Given the above arguments, the 
spinning NS is best approximated as a Jacobi-likc ellip- 
soid that is roughly rigidly rotating and, as we will show 
below, is actually very close to just being spheroidal. 

It should be noted that o ur model is decidedly diff er- 
ent than that employed by ICorsi fc Meszaroa (|2009( ) to 
follow NS spin evolution following a magnetar-powered 
long gamma-ray burst. They focus on spindown by grav- 
itational waves and dipole radiation and include no fall- 
back accretion, nor any magnetic or viscous effects on 
the internal fluid motions of the star. Because of this, 
their spin evolution asymptotes to a highly differentially 
rotating Dedekind ellipsoid instead. 

In this saturated state the NS emits GWs with an as- 
sociated energy loss of 



■ _ 32G^ 2 
gw ~ 5^~ { 1 ~ ' 



(12) 



where la = K n AIaf/5 are the components of the star's 
quadrupole, and a\ and 02 are the main axis of the tri- 
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axial figure. Angular momentum is removed at a rate 

A gw = E gw /n. (13) 

The corr esponding GW stra in amplitude as measured on 
Earth is ([Brady et alJfl998l ) 



2Gil 2 

"o — 4 n [hi - I22), 



(14) 



where D is the distance to the source. 

Given the above discussion, we summarize our proce- 
dure for estimating the GW emission from fallback ac- 
cretion as follows. 

(1) Accretion increases /3 on a timescale ~ 10 — 50 s un- 
til at a value /3 ps j3 c , GW production saturates the spin 
parameter. If the NS is destabilized due to gravitational 
radiation reaction, then f3 c w 0.16 — 0.18. But given the 
potential for viscous or accretion driven destabilization, 
we consider it reasonable to set f3 c ps /3 sec . 

(2) This /3 C implies an ellipticity e, spin frequency f2, 
and equatorial radius R e via equations ([?])■ l|8"j). and fj 10[) 
respectively. 

(3) For GW emission to maintain ps /3 C , it must re- 
move angular momentum at a rate that roughly balances 



the accretion torque, so that N t 



gw 



0. 



(4) Using this equilibrium condition, we can solve for 
In — I22 using equations (|12[) and (|13p . which can then 
be substituted into equation (fT4"| to find /irj. 

2.4. Estimates and Time- Dependent Evolutions 

To better demonstrate how our GW calculations will 
proceed, it is useful to present an example implemen- 
tation. Consider accretion onto an Mq = 1.3M@, 
Rq = 20 km NS. Setting f3 w /3 C and using /3 C = /3 soc , 
implies that the ellipticity and spin in the saturated 
state are e = 0.813 and n/2n = 542 Hz (n = 1) or 
0/27T = 457 Hz (n = 0.5). For n = 1, the mean 
radius is R = 1.09i?o and in turn the equatorial ra- 
dius is R e = 1.3LRo. For n = 0.5, i? = 1.035-R and 
R e = 1.24i?o- These values match what is plotted in 
Figure [2] at /3 = f3 scc (along the vertical dotted line), 
and are consistent with the general intuition that a more 
centrally concentrated mass distribution (i.e., higher n) 
is more strongly effected by rotation at a given f3. 

In Figure [3] we present example time evolutions of ac- 
creting NSs in the saturated state. The top two panels 
show the accretion rate from equation ([3]) using 77 = 1 
and the total mass given an initial mass of Mq = 1.3M0 
and using equation |QJ. It is assumed that the NS ac- 
cretes for 40 s before (3 = /3 C rit, and thus earlier times 
are not shown in this plot. We could potentially attempt 
to set the timescale at which the saturated state begins 
self-consistently with the early time accretion rate. The 
problem is that when this occurs will depend on the ini- 
tial spin of the core, which is difficult to estimate. We 
therefore leave the initial time for saturation as a free 
parameter and consider any time from ~ 10 — 50 s to 
be reasonable. As described previously, a given mass 
and value of /3 C is sufficient to solve for the spin rate (or 
GW frequency), which is plotted in the third panel of 
Figure [3J This depends on the density distribution of the 
star, and thus we plot / for both n = 1 (solid lines) and 
n = (dashed lines) polytropes. 
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Fig. 3. — Example frequency, mass, and strain amplitude evo- 
lution using the scheme described in these notes. The NS has 
Mo = 1.3Mq and R = 20 km, with both n = 1 (solid curves) and 
n = 0.5 (dashed curves), and rj = 1.0 for setting M. It is assumed 
that the NS accretes for 40 s during which < /3 cr it and no angu- 
lar momentum is lost. We then set f3 = /3 cr jt = 0.14 and follow 
the evolution up until M at 2.5Mq, at which point the NS would 
collapse to become a BH. 

When the GW torque is in equilibrium with the accre- 
tion torque {N gw + 7V acc ps 0), the quadrupole moment 
Q = hi - I22 is 



Q = 



5 (GMR e ) 1/2 Mc 5 



-1 1/2 



32 



Gtl 5 



(15) 



For M = 10 Mq s , we estimate a poloidal elliptic- 
ity e = Q/MRq = 7.0 x 10~ 3 [n = 1) or e = 1.1 x 10~ 2 
(n = 0.5). This demonstrates that the NSs we are consid- 
ering here are only deviating from sphericity by a small 
amount. This gives us some confidence that accretion 
will proceed unhindered by the triaxial shape of the NS. 
(In contrast, it is unlikely that accretion can occur onto 
a NS that is spinning fast enough for dynamical insta- 
bilities and has an extreme bar-mode shape.) It should 
be noted though that Q/MR 2 > 10 -5 exceeds the max- 
imum strain that a neutron crust can support (jShahaml 
1977). This crustal physics is not included in our calcu- 
lation, and a future, more detailed investigation should 
include the time-dependent effects of crust formation. 

Substituting this Q into the strain in equation (|14l) . 
and using a distance of D = 10 Mpc, gives ho = 
4.4 x 10~ 24 (n = 1) and 4.8 x 10" 24 (n = 0.5). We also 
plot the time-dependent strain amplitude in the bottom 
panel of Figure [3] This is similar, but slightly larger than 
these estimates since the NS and its spin frequency are 
increasing with time even as j3 remains fixed. 

3. DETECTABILITY OF SIMULATED WAVEFORMS 

In this section we survey upcoming GW interferome- 
ters and assess the detectability of GWs from a NS sub- 
ject to fallback accretion. Focusing on realistic excess 
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cross-power searches, we show that second-generation de- 
tectors such as Advanced LIGO (aLIGO) may be able 
probe accretion-powered NS models out to astrophysi- 
cally interesting distances: 17Mpc. 

3.1. Gravitational-wave interferometers 

Recent years have seen the development of a worldwide 
networ k of GW interferome ters. The initial LIGO exper- 
iment (|Abbott et al . 2009b) achieved design strain sensi- 
tivity: w2x 10~ 23 Hz 1 / 2 in the most sensitive frequency 
range of 100-200 Hz. The next generation of interfer- 
omet ers such as Advanced LIGO (aLIGO) ( Harry et al.l 
120101) . A dvanced Virgo (lAcernese et all (20061) . and 
KAGRA (jKuroda et al]|2010D are expected to begin tak- 
ing data starting in ss 2015, while the collection of 
data from GEO (jAbadie et al.1 120111 ) is ongoing. The 
upgraded aLIGO/ Virgo experiments are expected to 
achieve a factor of ten higher better strain sensitivity 
than initial LIGO/Virgo. In our simulation study below 
we show that aLIGO can probe fallback accretion signals 
to objects as distant as the Virgo Cluster. 

3.2. Waveform simulation 

In order to assess the detectability of fallback accretion 
events, we inject simulated waveforms into Monte Carlo 
noise and determine the distance at which signals can be 
distinguished from noise. Far from the source, we can 
write the metric perturbation as 



h + (t) = h (l + cos 2 i) cos $(t) 



and 



hx{t) = 2/ioCOS6sin$(t) 
where ho is calculated from equation (|T4|) . 



$(t) = / 2irf(t')dt', 
Jo 



(16) 



(17) 



(18) 



is the time-dependent phase, and t is the inclination an- 
gle, (which we later set to zero by assuming that the 
source is face onJJ We calculate h+(t), h x (t) by numer- 
ically evaluating the integral in equation 1181 

3.3. Injection recovery 

We envision an excess cross-power search (see 
IThrane et aHl20lU|P7estegard et al.ll2012[ ) for long-lived 
GWs associated with a well-localized electromagnetic 
counterpart, so that the direction and time of the signal 
are constrained, (though more computationally expen- 
sive "all-sky" searches are possible as well). We consider 
a simple network consisting of two 4 km initial LIGO ob- 
servatories operating at design sensitivity; one in Han- 
ford (HI) and one in Livingston (LI). 

We assume that the GW frequency of the signal is be- 
tween 700-2400 Hz, which is true for the models consid- 
ered here (see Table [1]). We assume that the electro- 
magnetic trigger constrains the on-source region (during 
which the signal may be present) to 1000 s in duration. 
Such an on source region is similar to the timescale over 



which the start of a supernova can be i dentified from 
shock breakout of a compact progenitor dNakar fc Saril 
I2010T) , and considerably longer than the typical timescale 
of a long gamma-ray burst. On the other hand, if the 
supernova progenitor is more extended (like a red super- 
giant), Hiemttie_c^i^oj^e£egion could be ~ 3 — 10 times 
larger (|Schawinski et al1l2008|) . We explore how sensitiv- 
ity varies with the on-source duration below. We assume 
that electromagnetic measurements constrain the direc- 
tion of the source to better than 0.17°, inside which at 
least 90% of the Hl-Ll point-spread function is contained 
for a source with / < 2400 Hz. (If the source cannot be 
this well localized, the sky can be tiled in order to search 
over multiple directions.) 

With in this 1000 s x 1700 Hz on-source region, and fol- 
lowing IThrane et all (|2011| ). we create a spectrogram 
of signal-to- noise ratio SNR(i;/), which is proportional 
to the cross-correlation of the HI and LI strain data. 
Here / refers to the frequency bin in a discrete Fourier 
transform centered on time t. We use 1 s-long, 50%- 
overlapping, Hann-windowcd data segments, which yield 
spectrograms with a resolution of 0.5 s x 1 Hz. Note that 



(SNR(t; /)) = hi (1 - cos 2 l) 



(19) 



3 Following [Prix fc K risiman j2009h . we work in a right-handed 
basis {ft, v, —h} where n points from the detector to the source, 
and fi, v are aligned the principal polarization axes of the GW. 



and that noise fluctuations can be both positive and neg- 
at ive. A forma l deriv ation of SNR(i; /) can be found 
in IThrane et al.l (|2011l ). An example of a SNR(i; /) spec- 
trogram can be seen in the two left panels of Figure [4] 

The problem of detection is to identify a track of 
positive- valued pixels in the p resence of noise. We em - 
ploy a track-search algorithm (|Prestegard et al.ll2012bl) . 
which looks for clusters of positive pixels^. The pixels 
are combined to determine the SNR for the entire clus- 
ter denoted SNR(c). Note that SNR(c) is distinct from 
SNR(i;/), which is the SNR associated with individual 
pixels in a spectrogram. We apply the track-search algo- 
rithm to Monte Carlo Gaussian noise in order to deter- 
mine the threshold for an event with false alarm proba- 
bility 0.1%. We find this threshold to be SNR(c) = 23. 

Real interferometer strain noise is only approximately 
Gaussian and includes noise bursts (glitches) and other 
non-stationary effects. As a rule of thumb, glitches 
have a tendency to push the required detection thresh- 
old higher than it would be in Gaussian n oise. How- 
ever, it was shown in lPrestegard et ail (|2012f) that initial 
LIGO cross-correlated data can be cleaned to the point 
where it is essentially indistinguishable from ideal Gaus- 
sian noise in the context of a search for long 0(> 10 s) 
narrowband GW transients. This is especially true in the 
shot-noise regime at high frequencies / > 500 Hz where 
our signal resides. Thus, we expect that the threshold 
of SNR(c) = 23 is representative of what can be achieved 
with realistic interferometer data. 

Next, we inject simulated signals into the Monte 
Carlo noise in order to determine the amplitude above 
which 50% of the signals are recovered above threshold. 
The signals are injected assuming an optimally oriented 
source (t = 0) and in the direction that maximizes the in- 
terferometer network sensitivity. (The sky-averaged sen- 
sitivity is « 2.5 less than the optimal sensitivity.) Plots 
illustrating the recovery of injected signals are presented 

4 For othe r examples of clustering algor ithms used in GW data 
analysis, see lKhan L Chatteriil |200g) and lRaffai et al.l l[2007ft . 
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TABLE 1 

Fallback accretion models with associated aLIGO Detection Distances' 1 



M max (M Q ) R (km) 77 / min -/ ma x (Hz) t dur (s) D bm . st b (Mpc) D m! c (Mpc) 



2.5 


20 


1 


1098-1512 


208 


9.5 


125 


2.9 


20 


1 


1098-1655 


354 


10.0 


135 


2.5 


15 


1 


1690-2328 


208 


5.0 


31 


2.5 


25 


1 


785-1082 


208 


17.0 


195 


2.5 


20 


10 


1174-1528 


37 


12.5 


115 


2.5 


20 


0.5 


1093-1508 


385 


8.5 


120 


2.5 


20 


0.3 


1091-1510 


641 


7.5 


145 


2.5 


20 


0.1 


1089-1510 


3100 


< 7.5 d 





a Distances are calculated for a false alarm probability < 0.1% and a fal se dismissal prob ability 
< 50% using Monte Carlo aLIGO noise (high-power, zero-detuning IIThe LSU 120101 )'). We 
assume an initial NS mass of 1.3Mq. Note, for comparison, that the Virgo Cluster is ~ 
16.5 Mpc away from Earth, while the Andromeda Galaxy is « 0.8 Mpc away. Sources are 
assumed to be optimally oriented. 

b Detection distance for a realistic excess cross-power search. 

c Detection distance for a highly idealized matched-filtering search. 

d Due to its length, more work beyond our present scope is required to carefully study this 
3100 s signal. We note that the detection distance is likely to be less than a signal with 
identical parameters except for a larger r\. When r\ is decreased, the power is spread out more 
diffusely over a longer duration, making the signal harder to detect. 




Fig. 4. — Recovered signals injected into aLIGO Monte Carlo noise: M = 1.3M Q , R = 20km, r) = 1, D = 5 Mpc (top) and M = 1.3M Q , 
R = 25km, T] = 1, D = 10 Mpc (bottom). Left: a spectrograms of SNR(i; /). We have cropped the spectrogram from its full size so that 
the signal is visible by eye. Right: candidate clusters identified by the track-search algorithm. (The small green and blue clusters in the 
top-right plot are noise fluctuations, while the long, thin, red track is the reconstructed GW signal.) The cluster SNR, SNR(c) = 168, 111 
(top, bottom), is well above the 0.1% false alarm probability threshold of 23. 
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in the two right panels of Figure SJ By varying the dis- 
tance to the source, we can determine the distance at 
which we can observe a signal with false alarm proba- 
bility < 0.1% and false dismissal probability 50%. This 
"detection distance" denoted -Dburst is given in Table [1] 
for different combinations of parameters. We also in- 
clude the corresponding distance that can be obtained 
with a highly idealized matched filtering analysis D m {. 
We return to the topic of matched filtering below. 

We find that, depending on the source parame- 
ters, aLIGO can observe signals out to distances of 
« 5 -17 Mpc. Thus, aLIGO will be able to probe 
fallback-powered GW emission from NSs in nearby galax- 
ies such as Andromeda (0.8 Mpc) and perhaps as far away 
as the Virgo Cluster (16.5 Mpc). These distance esti- 
mates are obtained using an existing pattern recognition 
search algorithm. It is possible that the sensitivity can 
be improved by developing a specially tailored search al- 
gorithm. 

Until now we have assumed that the GW signal can 
be confined to a 1000 s on-source region. In order to 
determine how this assumption affects the search sensi- 
tivity, we also consider a longer 3000 s on-source region. 
We find that the 0.1% false alarm probability thresh- 
old grows from SNR(c) = 23 to SNR = 25. Repeat- 
ing the above injection study, we find that the detection 
distance for a 3000 s window is about 3% smaller than 
for a 1000 s window — e.g., the detection distance for the 
R = 25 km waveform shown in Figure |4] changes from 
17.0 Mpc to 16.5 Mpc. Thus, the falloff in sensitivity is 
modest as the on-source region is increased in size, which 
in turn suggests that the analysis we have sketched out 
can be performed for a wide range of electromagnetic 
triggers. In the next section, we investigate the best 
theoretically possible (though probably unrealistic) sen- 
sitivity that can be achieved with matched filtering. 

3.4. Matched filtering estimates 

For GW signals of a known form, the theoretically 
optimal searc h strategy is matched filtering (see, e.g ., 
lAbadie et afl l2012al |bT lOwen fc Sathvaprakashl Il999f ). 
The fallback accretion powered signals we consider here, 
however, will prove challenging to target with a matched 
filtering search. Such a search must include a template 
bank that spans the space of possible signals. An incom- 
plete template bank can result in faulty upper limits if 
the true signal falls outside the template space. In order 
to generate a complete template bank, we require firm 
knowledge of the details of the waveform's phase evolu- 
tion. In presenting the above model, we have not aspired 
to this degree of accuracy. Moreover, even if a complete 
and accurate model could be written down, there may be 
computational challenges associated with performing the 
search, especially for long signals with many parameters. 

Nonetheless, it is useful to compare the detection dis- 
tances calculated using the excess-power technique to es- 
timates for what can be achieved with matched filtering 
as this places an upper limit on the detection distance 
that can be achieved through improvements to the data- 
analysis scheme. Rather than create a bank of templates 
to span the space of our signal, we assume that every 
signal parameter is known including: the start time of 
the signal, the initial phase, the initial mass Mo, the ac- 
cretion parameter rj, the polytrope index, the direction, 



and the critical value of (3. As before, we assume optimal 
orientation of the source, optimal orientation of the de- 
tector network, a false-alarm probability of 0.1% and a 
false dismissal probability of 50%. This highly idealized 
analysis allows us to use just a single template. We find 
that highly idealized matched filtering allows us to ex- 
tend the detection distance by a factor of 6 to 19. This 
factor varies depending in part on the efficiency of the 
pattern recognition algorithm for different signal types. 

3.5. Other Gravitational- Wave Signatures 

Besides the GWs from fallback, there may be other 
GW signals associated with supernovae and the creation 
of a BH. A variety of mechanisms have been proposed 
for GW emission from supernovae, including the acous- 
tic mechanism, rotational instabilities, rotating collapse 
and bounce, c onvection a nd standing accretion shock in- 
stabilities (see lOttl 120091 and references therein). Com- 
peting models predict GW fluence that spans six orders 
of magnitude (2 to 2 x 10 6 erg cm -2 for a galactic source 
at lkpc), but at least one model (the acoustic mecha- 
nism) predicts strain amp litude high enough for aLIGO 
to observe in Andromeda (|Leonorj|20l"oHOtt et al.ll2006t) . 
It is therefore possible that that the fallback GWs we 
describe here could potentially be present from the same 
event that produces a short GW burst in Andromeda. 

In addition to the precursor signal from core collapse, 
it is interesting to consider the possibility of a postcur- 
sor signal from the ringing black hole eventually forme d 
through accretion (|Frver et all 120021: lOtt et al.l [20111) . 
There are significant uncertainties surrounding the strain 
amplitude of the black hole ringdown signal following 
accretion-induced collapse, b ut the signal mor phology 
is relatively well constrained (jFrver et al.ll2002[ ). For a 
newly formed black hole with mass 2.5M Q , we expect a 
damped sinusoid sign al with frequency f ~ 6.5 kHz and 
quality factor Q « 3 (| Abbott et al.ll2009a|) . The strain 
amplitude depends on the fraction e of energy radiated in 
GWs and the mass /x falling into t he black hole, the reby 
exciting the quasinormal modes ()Frver et al.ll2002T) . If 
we optimistically assume e = 0.01 and [i = 0.1M Q , the 
strain amplitude is ho w 6 x 10 -22 , which implies a 
matched filtering detection distance of 0.15 Mpc (assum- 
ing an optimally oriented source, requiring a false alarm 
probability = 0.1%, and a false dismissal probability = 
50%). Thus the possibility of observing a ringdown sig- 
nal is probably limited to sources in our own galaxy. 

4. ELECTROMAGNETIC COUNTERPARTS AND RATES 

Given that fallback accretion onto a NS can be ob- 
served with aLIGO out to a distance of ~ 17 Mpc, it 
is worth exploring the broader context of what kind of 
events could be associated with potential GW detections. 
It has been well-established in the case of compact merg- 
ers that joint GW and electromagnetic detections pro- 
vide the greatest opportunity for learning about these 
systems. This should likewise be true for core-collapse 
supernovae. In the following discussion we draw upon 
recent theoretical and observational progress in under- 
standing core-collapse supernovae to address potential 
electromagnetic counterparts and estimate rates. 

To get fallback accretion as needed for GW produc- 
tion requires a few key ingredients: (1) sufficient angular 
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momentum for disk formation, (2) a sufficiently low ex- 
plosion energy to not unbind the entire envelope, and 
(3) a progenitor that is sufficiently massive and compact 
for eventual BH formation. Disk formation within col- 
lapsing stars has traditionally been the focus of studies 
trying to produce long gamm a-ray bursts. For exam- 
ple, a large part of the work by iMacFadven et all (|2001h . 
which motivated the fallback model we employed, was 
presented to study the hyperaccretion that would occur 
following BH formation. Since we too require rotation to 
spin up the NS and generate GWs, this may mean that 
hyperaccretion and jet formation may follow. But this is 
not certain, because the exact conditions required for jet 
formation are not well-known. We therefore divide our 
discussion between events that have jets and those that 
do not. 

4.1. Counterparts of Events with Jet Formation 

Even if a jet is produced, the observational signature 
can vary greatly depending on whether the jet breaks out 
of the envelope. For this to occur, a substantial fraction 
of the envelope must be removed. This is consistent with 
the correlat ion of long gamma-ray b ursts with Type Ic 
supernovae (jWooslev fe Bloo m 2006), which do not have 
signatures of hydrogen or helium. If there is a sufficiently 
large hydrogen envelope to prevent the jet from cleanly 
breaking out, a jct-powcred Type II supernovae may re- 
sult instead. Such an event might look like SN 2010jp 
(jSmith et al.ll2012l ). which showed evidence for a bipolar 
outflow. In addition, it had a low output of radioactive 
56 Ni (< 0.003Mq), consistent with appreciable fallback, 
and signs of low meta llicity, consistent with less mass 
loss (jHeger et al.ll2003l) . A jet-driven explosion was also 
suggested for SN 200 1ig based on spectropolarimetry and 
line-p rofile evidence (jMaund et al.ll2007tlSilverman et al.l 
I2009T ) . Other studies have attempted to use polariza- 
tion to infer the amount of asymmetry in the supernova 
(whi ch could be driven by a jet) , but this is less conclu- 
sive (|Leonard et al.l l2000l 120011) . It is enticing though 
that the amount of asymmetry seems to be inversely 
relate d to the hydrogen envelope mass (jChornock et al.1 
1201 ll and references therein), which is consistent with 
the generally expected picture. If the hydrogen envelope 
buries all signs of a jet, even a weak amount of subse- 
quent hydrogen- reco mbination may produc e a luminous 
red nova-like event dKulkarni et all 120071 : iBond et all 
20091: ITnompson et al.ll2009h . as has been speculated by 



TABLE 2 

Recent Nearby Core-collapse Supernovae 51 
Name Type Galaxy Distance (Mpc) 



Quataert fc Kasenl (j2012l ) 



4.2. Events without Jet Formation 

If the majority of the material just falls into the BH 
or if there is little 56 Ni production, the electromagnetic 
signature may be ver y weak, producing an "unnova" 
(|Kochanek et al.ir2008|) . The prospect of identifying dy- 
ing stars that cannot be observed with traditional meth- 
ods would be an exciting use of GWs. On the other hand, 
there is currently little evidence for such events, and no 
apparent need of them to explain discrepancies between 
rates of different s upernova types ver sus a standard ini- 
tial mass functions (|Smith et al.ll2011[). Furtherm ore, the 
nearby Type II-L SN 1979C (jPanagia et al.lH980h clearly 
had a bright, explosive d isplay, but some have argued 
that it left a BH remnant ([Patnaude et al.ll2011l ). A dim 



2012aw 
2012A 
2011ja 
2011dh 
2009dd 
2008bk 
2008ax 
2007it 
2007gr 



II-P 
II-P 
II-P 
lib 

II 

II 

lb 

II 
Ib/c 



M95 
NGC 3239 
NGC 4945 

M51 
NGC 4088 
NGC 7793 
NGC 4490 
NGC 5530 
NGC 1058 



10 

8 

3.6 
7.1 
15.8 
3.9 
15 
NA b 
8.4 



a Summarized from lists of the Astronomy Sec- 
tion^ Rochester Academy of Science (ASRAS), 
http: / /www. rochcstcrastronomy.org/snimagcs/ 

T5 - 



The distance for this event has not been re- 
ported. Nevertheless, since it was the brightest 
supernova of 2007 (apparent magnitude), it was 
likely within the sa 17 Mpc that we require. 

or nonexistent event also seems to be disfavored on phys- 
ical grounds. The outer envelope of a red supergiant is 
very loosely bound, and just a weak shock will eject it 
and produce a sign ature from hydrogen recombination 
(jDessart et al.ll2010D . even if all the radioactive material 
falls into the BH. This would argue that a supernova 
with fallback may produce a Type II-P like light curve, 
even if it is relatively weak. Again, this might lead to a 
connection with some subset of the luminous red novae. 

4.3. Rate Estimates 

Since it has been difficult to definitively connect any 
electromagnetic transients with fallback accretion and 
BH formation, the rate of such events, especially within 
the « 17 Mpc distance needed for GW detection, is un- 
certain. Long gamma-ray bursts may be ideal sources, 
since GWs from fallback may be a precursor to jet for- 
mation. Such a detection would provide important ev- 
idence that some long gamma-ray bursts are powered 
by accretion versus m odels that favor young magnctars 
(|Metzger et al.ll2011l and references therein). Unfortu- 
nately, the closest recorde d long gamma-ray b urst was 
at a distance of 37 Mpc (|Galama et al.l [T998|) . making 
the prospect of a nearby long gamma-ray burst unli kely. 
Based on rate estimates from (jChapman et al.ll2007l ). we 
expect 3kyr _1 within 17 Mpc. 

On the other hand, ~ 1 — 2 supernovae typically oc- 
cur within 17 Mpc each year; see Tabled] Most of these 
nearby events are the common Type II-P supernovae, 
which are perhaps the least likely to have fallback accre- 
tion, although there are a few other types. Nevertheless, 
given the uncertainties in the electromagnetic signature 
from fallback, events like these would be worth checking 
with next generation gravitational wave interferometers. 
The Hubble Space Telescope can detect progenitors of 
these events i n pre-explosion imaging out to ps 30 Mpc 
()Smarttl 12009). So in principle, we might expect that 
any supernova that has GWs detectable from fallback 
will therefore have an identified progenitor. In practice, 
only ~ 1/3 of the time is the archival imaging available 
for the location of the nearby supernovae. Also problem- 
atic is that there still has never been a progenitor star 
identified for a Type lb or Ic supernovae, which remains 
an outstanding puzzle. 

Although much more speculative, if there is fallback 
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accretion associated with some subset of luminous red 
novae, this could provide an important increase in the 
rate of fallback GW detections by ~ 50% or more 
since these events are relatively nearby ([Thompson et ahl 
Since luminous red novae are often shrouded in 
dust, they can be difficult to study and may even be 
missed electromagnetically in some instances. GWs pro- 
vide a method to potentially detect and learn about these 
mysterious transients. 

The rate of events without any electromagnetic signa- 
ture is the most difficult to estimate. Monitoring data 
for M81 from the Large Binocular Telescope formally set 
the rate as < 80 times tha t of normal SNe at 90% con- 
fidence (jPrieto et alj|2008t) . which is not especially lim- 
iting. Perhaps more constraining, the similarity of mas- 
sive star formation rates and supernova rates and the 
nondctcction of a diffuse supernova neutrino background 
both indicate that the rate of events without electromag- 
netic signatures cannot significantl y exceed the rate of 
obser ved core-collapse supernovae (jHopkins k, Beaco~ml 
120061) . Even though a systematic surve y is necessary to 
provi de stringent limits on such events (jKochanek et al.l 
12008( 1. the prospect that GWs could help find a missing 
class of dying stars is worth exploring in more detail. 

5. DISCUSSION AND CONCLUSION 

We investigated the production of GWs in a subset 
of supernovae where a NS is produced in the core, but 
which subsequently collapses to a BH after pa 30 — 3000 s 
of fallback accretion. Using semi-analytic models for the 
spin evolution of the NS, we estimate typical GW fre- 
quencies of ps 700 — 2400 Hz as well as the corresponding 
strain amplitudes. We estimate that such events will be 
detectable by Advanced LIGO out to pa 17Mpc. Given 
the messy process of fallback accretion, we argue that 
this figure is much more robust than estimates made with 
an overly optimistic matched-filtering technique. Finally, 
we discuss a number of electromagnetic events that may 
be associated with this unique GW signal. Although 
such connections are still very uncertain, there are likely 
^ 1 — 2 supernovae events per year that will be worth 
searching for fallback GWs. 



In this initial work a number of simplifying assump- 
tions were made in the GW model. An important focus 
of future studies is to investigate more realistic progen- 
itors. Future improvements may include: using realistic 
fallback rates from numerical models, resolving the time- 
dependent accretion disk that mediates the flow of mat- 
ter onto the NS, and better understanding what value of 
j3 c is appropriate for a NS being spun up by fallback. 
Furthermore, a better survey of the range of fallback 
rates and expected timescales should be conducted, es- 
pecially in the case of stellar models that may have been 
ignored because they do not produce long gamma-ray 
bursts. Binary interactions may also play an important 
role in generating the needed rotation rates (for exam- 
ple, see lWooslev fc Hegerll2012l ). a fact that is maybe not 
surprising given that the majority of ma ssive stars are in 
close binary systems (jChini et alJl2012j ). 

In the event of a detection, there are a number of in- 
teresting things we might learn from the GW signature. 
For a given frequency, there is a one-to-one correspon- 
dence with the density of the NS, given constraints 
on the polytropic index n and /3 C via equation (j8[). 
Therefore determining the frequency at the moment 
before collapse to a BH would directly constrain the 
density at the maximum mass of a NS, an important 
const raint on its equation of state (|Lattimer fc Prakashl 
2007). The main uncertainty is in j3 c . We have ar- 
gued that f3 c « f3 scc is the most likely possibility, but 
even if f3 c pa 0.1 — 0.15, this would only introduce a 
pa 30% error in the density assuming the frequency is 
precisely measured. Nevertheless, this highlights the 
importance of determining /3 C in future, theoretical work. 

We thank Christian Ott for detailed comments. We 
also thank Lee Lindblom for helpful discussions on non- 
axisymmctric instabilities of rotating NSs, and Andrew 
MacFadyen for feedback on the observational impact of 
fallback accretion during supernovae. A.L.P. was sup- 
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